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ABSTRACT 


The study outlines a procedure for adapting existing hyper- 
sonic flow field programs to nonequilibrium conditions. The 
procedure treats bound-bound transitions between electronic 
levels in molecular species in which electron impact is 
considered to be the dominant excitation mechanism. The 
formulation of expressions for excited state population and the 
radiative transfer equation are presented. The formulation 
indicates how for nonequilibrium of the radiative transfer 
equation can be written in terms of spectral source functions 
and absorption coefficients so that existing programs for 
radiative properties can be used as well as existing radiative 
transfer solution programs . 
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-Einstein coefficient for spontaneous emission 2 -* 1 
-Einstein coefficient for stimulated emission 2 -» 1 
-Einstein coefficient for stimulated absorption 1 -*• 2 
-Planck function at frequency u and temperature T 
-velocity of light 

-fraction of the particles in an electronic state which are 
in the i'th vibrational state 
-statistical weight factor of state i'th 
-Planck constant 

-specific intensity at frequency v 

-mean intensity at v specific intensity averaged over 
direction 

-Boltzmann constant 

-concentration of particles in lower electronic state - 
number/unit volume 

-concentration of particles in upper electronic state - 
number/unit volume 

-concentration of electrons - number/unit volume 
-vibrational partition function for the upper electronic 
state evaluated at temperature T 
-vibrational partition function for the lower electronic 
state evaluated at temperature T 
-Franck-Condon factor for a transition between the v' ' 
vibrational level of the lower electronic state and. the v' 
vibrational level of the upper electronic state 
-Source function at frequency v 
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-translational temperature of free electrons 
-Boltzmann temperature of the lower electronic state 
-Boltzmann temperature of the upper electronic state 
-Boltzmann vibrational temperature of the lower electronic 
state 

-Boltzmann vibrational temperature of the upper electronic 
state 

-vibrational quantum number of the upper electronic state 
-vibrational quantum number of the lower electronic state 
-depth coordinate 

-energy of the v' vibrational state relative to v' =0 
-energy of the v' ' vibrational state relative to v" =0 
-energy level of the upper electronic state v' =0 
-energy level of the lower electronic state v' ' = 0 
-energy difference between electronic states - e-^) 

-cosine of polar angle between intensity and depth 
coordinate 

-absorption coefficient per unit depth at frequency u 
-frequency 

-optical thickness at frequency u 

-normalized line shape profile for lower electronic state 
-normalized line shape profile for upper electronic state 
-collisional excitation rate 1 -*• 2 
-collisional quenching rate 2 -+ 1 

-identifies an i'th electronic transition and/or specie 
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INTRODUCTION 


Nonequilibrium occurs in a gas when an imbalance in excitation and 
quenching (de-excitation) processes results in state populations which 
deviate from the Boltzmann distribution characterized by the gas 
temperature. In the general case, the excitation and quenching processes 
include both kinetic (collisional) and, radiative processes. For the 
molecular species, the subject of this study, each state is specified by 
a rotational quantum number, a vibrational quantum number, and the 
electronic level. Under the most sever nonequilibrium condition, each 
state is characterized by a different temperature which is not the 
translational temperature of the gas; in fact, the translational 
temperatures of the species are not necessarily equal. 

A severe perturbation of atmospheric gases by an entry shock results 
in a nonequilibrium layer of gas adjacent to the shock. In this 
nonequilibrium layer, the rapid rise in random translational energy of 
the post-shock gas is re-distributed by kinetic processes into the 
internal states of the gas (rotational, vibrational, electronic, chemical 
change , and ionization) . 

The combination of hypersonic entry velocities with the low densities 
encountered in upper atmospheric entry produces a shock-layer 
characterized by high temperatures and low collision rates. For these 
conditions, the nonequilibrium region is extensive and the prediction of 
radiative transfer from this region is essential in prediction of the 
instantaneous heating from the shock-layer. 

Nonequilibrium radiative heating is encountered by every body during 
the early portions of entry. During this portion of the trajectory, 
radiative heating of the body poses a more severe and uncertain problem 
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than conventional equilibrium heating estimates would indicate. The 
problem was studied in the Apollo program. For Apollo, the 
nonequilibrium heating did not pose a severe problem in total system 
thermal management as Apollo plunged quickly through those conditions for 
which nonequilibrium dominates. A re-examination of the early flight and 
laboratory data has been presented by Sutton [1]. 

To date, the hypersonic shock layer heating has been handled by using 
either a chemical equilibrium or a chemical nonequilibrium flow field 
analysis to describe the spatial distribution of radiating species. An 
example of such a program (NEQ-VIS) is the viscous-shock-layer program of 
Moss [2] which has been extended for use in chemical nonequilibrium flow 
field analysis [3] . These applications are for conditions in which the 
collision rates are sufficiently high to assure kinetic equilibrium of 
the internal states to a common temperature and to assure kinetic rates 
exceed those of the radiative processes. For these conditions, the 
prediction of radiative transfer through the shock-layer follows the same 
procedure whether chemical nonequilibrium or equilibrium is considered, 
and programs such as Nicolet's (RAD) [4] can be used with the spatial 
distribution of temperature and radiating specie from the flow field 
programs . 

Systems of current interest [5], Orbital Transfer Vehicles (OTV) or 
(AOTV) , are proposed to operate in an entry-exit maneuver - "skipping" 
through the rarefied upper atmosphere. As a result, an OTV would dwell 
for a significant portion of its entry-exit traj ectory 'in conditions for 
which nonequilibrium dominates [6], As an example, consider an OTV entry 
of 10 to 11 km/sec which attains a minimum altitude of 80 km. For these 
conditions, the number density in the shock layer will be of the order of 
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10 particles/cm (number density at STP is 10 particles/cm ) and the 
collision frequency would be 10^ collisions/sec. Clearly, the collision 
quenching rate can be no larger than this; however, the radiative 
quenching rate can be several orders of magnitude larger. Thus, both 
radiative and kinetic rates govern state populations and nonequilibrium 
predictions of state populations are required. 

The situation where kinetic rates are dominant but do not balance can 
be termed kinetic nonequilibrium. For such situations, one solves the 
kinetic rate equations to determine state populations. With the 
populations, one can solve the radiative transfer problem using 
procedures similar to those for the equilibrium. By contrast, those 
situations associate with OTV's in which radiative rates are significant 
the radiative terms must be considered in the rate equations. 

To answer the question of the importance of nonequilibrium radiative 
transfer, one should solve the system of coupled rate equations for the 
state populations together with momentum, energy and radiative transfer 
equations for the complex shock-layer geometry of the OTV. This is a 
formidable problem. One approach, taken by Park [7], is to meticulously 
assess the rate processes in the context of a simplified flow field 
geometry. 

Below, we approach the problem with the question: Can programs under 
development or existing programs similar to NEQ-VIZ and RAD be modified 
to solve the OTV problem to conserve the considerable effort required in 
code of development? The logic behind the question is that the form of 
the continuum conservation equation for mass , momentum and energy and 
their boundary conditions are the same regardless of the state of 
nonequilibrium. The treatment of chemical rate processes would seem 
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extendible to kinetic rate processes. Similarly, if the nonequilibrium 
radiative transfer problem can be stated in terms of spectral absorption 
coefficients and source functions, then the spectral and spatial 
integration procedure of existing computer codes can be used. 

To explore the limitations associated with this approach, the 
character of a generic flow-field and radiative codes, will be briefly 
reviewed in the next section. A simple two level system is used to 
illustrate the conditions of nonequilibrium for the state population rate 
equation and the radiative transfer equation and to examine 
simplifications. The simplifications are then incorporated into general 
multi-state rate and radiative transfer equations. The result of this 
study is a procedure which enables one to modify existing codes for 
application to the nonequilibrium shock-layer of OTV's. 
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REVIEW OF CODES 


The computational problem involving an entry shock-layer can be 
viewed in three, not necessarily equal parts: 

Part (1) Satisfying conservation equations, boundary conditions and 
chemistry 

Part (2) Computation of state populations for radiative properties 

Part (3) Satisfy the radiative transfer equation 

Usually, these parts are considered sequential; however, with 
radiative terms coupling all three together, a sequential solution must 
be considered as iterative. 

The governing equations for the continuum regime of an OTV trajectory 
have been reviewed recently be Lee [8]. Comparing these equations, using 
order of magnitude arguments, with those incorporated into a program such 
as NEQ-VIS , one finds the governing equations satisfied to better than 
first order. This is not a surprise as the form of the conservation of 
mass, momentum and energy equations are basically the same, equilibrium 
or nonequilibrium. Certainly, kinetic nonequilibrium perturbs gas 
properties such as viscosity, conductivity, and specific heat; and these 
can be corrected, if one desires, in iterative cycles. 

More important are three aspects of Part (1) which are not 
incorporated into the generic flow-field codes with nonequilibrium 
chemistry. These are: 

(1) Computation of the translational temperature, T , for the 
electrons 

(2) Computation of the vibrational temperatures, T , for the 
electronic states 
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(3) Proper use of heavy particle translation temperature T, electron 
temperature T , and vibrational temperature T in the evaluation 

i 

of the chemistry rate constants. 

These three aspects can best be handled by the addition of the energy 

equation for electrons and rate equations for the vibrational states to 

the set of simultaneous equations treated by the program which handles 

Part (1) . Adequate treatment of in NEQ-VIS has been reported by Shinn 

[9], The importance of incorporating vibrational temperatures is at this 

point, arguable. It can be done; however, in some of Park's work, 

arguments have been presented for T^ to T^ equilibrium. The proper use 

of T, T and T in the evaluation of chemical rate constants has been 
v e 

treated recently by Park [10]. Thus, Part (1) can be handled by existing 
programs with: 

(1) the addition of an energy equation for electrons 

(2) proper use of T g , T and T v in rate constant evaluations 

(3) proper evaluation of radiative contribution in the energy 
equation 

Part (1) will supply a flow field map of T , T, chemical specie 

concentration, and T 's. 

v 

The Part (2) computation of state populations is treated in a later 
section of this study. At this point, we should note that for the 
temperature and the level of ionization expected, the primary excitation 
mechanism will be from the ground electronic state and due to electron 
impact excitation. Due to their high speed, electrons are likely to 
dominate heavy particles in excitation. The excitations of rotational 
and vibrational states are sufficient fast to assure a Boltzmann 
distribution for each rotational and vibrational system. However, the 
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Boltzmann vibrational temperature for each electronic state may be 
different. Essentially, Part (2) converts the spatial distribution of 
T , chemical specie number density, and N g into a spatial distribution of 
electronic state populations fractions or temperatures. In conclusion, 
recall that the stimulated emission and absorption processes result in a 
coupling of Part (2) to the solution for Part (3) and that for the 
assumption of an optically thin gas these stimulated terms are neglected. 

For Part (3) , codes such as RAD handle radiative transfer as two 
problems. One consists of computing the optical properties (absorption 
coefficients and source functions) for every spatial point based upon the 
temperature and number density field supplied. The second problem 
involves the integration of the radiative transfer equation to achieve 
+he spatial and spectral values of the radiative intensity or of the 
integrals of the intensity. From these quantities, radiative heating of 
the OTV can be determined as well as the radiative coupling terms which 
should be used for the next iteration through Parts (1) and (2) . 

The procedure for solving the radiative transfer equation and 
determining heating rates is independent of nonequilibrium considerations 
providing the radiative transfer equation for nonequilibrium can be 
written in the standard form involving a spectral absorption coefficient 
and source function. The problem posed by the standard form of the 
radiative transfer equation is not so much in the absorption coefficient 
as in the source function. It is clear that each absorption coefficient 
will contribute to the frequency interval associated with its electronic 
state and it is clear how to sum the absorption coefficients spectrally. 
What is not clear is: What does one replace the Planck source function 
with to achieve the classic form of the radiative transfer equation? 
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This source function and its logical development is treated in a later 
section devoted to the radiative transfer equation. It will be shown 
that taking simple combinations of absorption coefficients and Planck 
functions evaluated at the nonequilibrium state temperatures will enable 
the spectral and spatial solution procedure of the RAD code to operate as 
if solving an equilibrium problem with spatial variation of absorption 
coefficients and source function. 

Returning to the former problem of optical properties, several 
approaches are available. First, a line by line computation has been 
used by several investigators . Such an approach requires a knowledge of 
the line broadening mechanisms to be used with an extensive data base and 
accounting procedures. Second, a JORL "just overlapping-rotational-line" 
band model [11] and [12] has been used extensively. This model portrays 
a simple variation in absorption coefficient within a vibrational band by 
summing over each rotational line which is assigned a line width 
conforming to the line spacing. A third procedure which is referred to 
as the "bandless model" or average absorption coefficient model considers 
the coefficient to be constant in and evaluated by averaging over a 
frequency interval Ai/. In this way, sums over all rotational lines lead 
to expressions which are independent of rotational temperature and 
parameters. This procedure was used in the early work of Meyerott [13], 
with Av = 6000 cm ^ and the frequently cited work of Biberman [14] with 
Ai / = 3000 cm \ These frequency intervals would sum together several 
vibrational transitions. To be consistent with the average coefficient 
assumption the spectral variation should be portrayed as a series of 
steps; however, the variation of amplitude is frequently represented by a 
smooth curve with the amplitude associated with the center of the 
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frequency interval. This model and procedure are incorporated into RAD 
[4] and will be used in a later section in which the none qui librium 
radiative transfer equation is examined. The result of this examination 
indicates that the nonequilibrium optical properties can be computed to 
the same accuracy as equilibrium optical properties using the procedures 
of RAD providing the properties are evaluated at the proper 
nonequilibrium temperatures supplied by Part (2). Thus, the use of 
proper sums and temperatures allows RAD to be used without major 
modifications . 
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TWO -LEVEL ILLUSTRATION 


A two-level system similar to that used in an earlier investigation 
[15], provides a simple illustration of the problem of predicting state 
populations for nonequilibrium and for observing the consequences of 
these populations on the radiative transfer problem. This will allow 
critical points to be observed uncluttered by the complex notation 
inherent in a multi-level system. Later sections will apply the findings 
of this section to more complex multi-level systems. 

For a radiating specie possessing two nondegenerate states between 
which both radiative and kinetic transitions occur, the rate equation and 
the radiative transfer equation take the form: 

^ =* -V A n + B 2i j + N n n + Vn > (1) 

+ n i( b 12 J + n o<$> +M e n<2 ) > 

dl 

M ~ = -<hi/*/4ir) [N 1 B 12 I i/ -N 2 (A 21 + B 21 I^) ] (2) 

dz 

The terms of the RHS of Eq. (1) are the rate of quenching of the 
population N 2 of the upper state and the rate of excitation from the 
lower state. The quenching processes are spontaneous emission, 
stimulated emission, collisional de-excitation by heavy particles, and 
collisional de-excitation by electrons. The excitation processes are 
stimulated absorption, excitation by heavy particles, and excitation by 
electrons. The differentiation is the material derivative associated 
with the stream lines of the flow. The differentiation of Eq. (2) is 
along a ray path. 
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Rearranging Eq . (2) yields the classical radiative transfer equation: 


dl 

M— = (hi/*/4ir) (N 1 B 12 - N 2 B 21 ) [S^ - IJ (2b) 

dz 

where the source function 

S, = N 2 A 21 /(N 1 B 12~ N 2 B 21 ) = (2h^ 3 /c 2 )/[ (Nj/N^-l] 

3 2 

using B^ 2 = B^ and A 21^21 = / c ’ anc * ° ne reco 6 n i- zes the 
absorption coefficient as: 

= (hi /*/4ir) [(N^d - N 2 /N x )] 

Now for many applications the optically thin assumption is used. 

This is based upon the assumption that the optical depth is small. 

For a uniform gas, the optical depth is proportional to the produce of 
the absorption coefficient and the ray path length. Recall that the 
absorption coefficient is inversely proportional to the photon mean free 
path. Thus, large mean free paths compared to the characteristic ray 
path dimensions implies small optical depths and photons escape the 
emitting volume without participation in the stimulated absorption or 
emission processes. 

For the optically thin assumption, Eq. (2b) is conventionally written 

as : 


dl 

v 

M 

dz 


= M. 


S 

u 


(3) 
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implying that the stimulated absorption and emission processes are - 
negligible. Reconsidering the definition of the absorption coefficient 
and source function yields a more useful form of Eq. (3) for the 
optically thin gas. 


dl 

H “ = (3b) 

dz 

The quantity hiM^^ is the radiative power emitted per unit volume 
and will be noted by Q/v. 

Returning to Eq. (1) for the population of the upper state ^ and 
seeking the steady state population i.e., the maximum population to which 
the process relaxes, we have set Eq. (1) equal to 0 to achieve the steady 
state result. 


N„ 


N, 


(e) 


N e°12 


+ B 12 J(1 - n 2 /n x ) 


N e°2t ) 


+ A 


21 


(4) 


where we have neglected the heavy particle collision terms based upon the 
argument that the columbic interaction by the electrons results in cross 
sections larger than heavy particle cross sections by an order of 
magnitude, that the speed of the electrons is over 250 times those of the 
heavy particles and that for a degree of ionization of 1% or more, the 
electron impact process exceeds the heavy particle process by more than 
an order of magnitude. Clearly, the heavy particle kinetics could have 
been retained; however, we would require the rate constants for them. 
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Now, if the numerator and denominator of Eq. (4) were divided by 
we would find the term [J (1 “ to be of the order of the 

terms which were considered negligible in the optically thin 
approximation. Thus, we shall neglect these stimulated processes and re- 
write Eq. (4) as: 



[ fl 12 N e 1 


°12 

N 2 

A 21 , 


. °21 


*■ 

n, „ 1 


0. J 'j 


12 

i 

12 e 

. 

°21 

~r 

. A 21 . 


< n i2V A 21> 

1 + (n 21 N e /A 21> 


Both ar *d ^21 are f unct *- ons t ^ ie temperature of the electrons and by 
the usual kinetic arguments relating forward and reverse rates: 

n 12 /n 21 = (N 2 /N l)eq - exp <-« 12 /W,> (6) 

Before examining Eq. (5) further, let us compare Eqs . (4) and (5) and 
ask the question: how does neglecting the stimulated absorption and 
emission terms influence the state population? Clearly, the dominant 
absorption process is an excitation mechanism which tends to increase the 
upper state population and so estimates in which these stimulated 
absorption terms are neglected, such as Eq. (5), will tend to under 
predict the excited state population. 

Returning to Eq . (5), notice that is the ratio of the 
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rate of kinetic quenching to the rate of radiative quenching and that at 
large electron concentrations, kinetic quenching will dominate and at low 
electron concentrations, radiative quenching will dominate. 

For large electron concentrations , i.e., Eq. (5) 

reduces to: 


(n !2 N e/ A 21 ) 
N /N = — ' — ~ 

V 1 (fl 21 N e /A 21 ) 


°12 

— = exp (-, 12 AT e ) 


= ‘Weq 


(7) 


using Eq. (6). This is the equilibrium solution which is maintained in 
equilibrium with the electrons by balancing of kinetic excitation and 
kinetic quenching. 

As the electron concentration is reduced, the kinetic quenching rate 
is reduced until it is below the radiative quenching rate. Thus, for 
N g « ^21^21 (5) reduces to: 


V N 1 = fl 12 N e /A 21 


( 8 ) 


where excitation by electron collision is balanced by radiative 
quenching. This solution is frequently referred to as "collision 
limiting . " 

Clearly for any prediction, Eq . (5) would be used with Eq. (6). 

However, it is useful to keep these two limiting solutions in mind. To 

achieve some understanding of the magnitudes involved in these two limits 

* 

consider the conditions for / A 2l = ^ 
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(N e ) " A 21 /n 21 ~ (A 21 /n i2 ) (n !2 /n 21 ) 


( 9 ) 


In Figure 1, the radiative lifetime and excitation rate data from [16] 

and [17] have been used to compute the "limiting electron concentration 

(N )* for 4 band systems: the first-positive (1+) , the second-positive 
e 

(2+) , and the Lyman-Birge-Hopf ield (LBH) of Nitrogen and the Schumann- 
Runge (SR) of Oxygen. The radiative lifetime is that of the optical 
transition while the excitation is always the resonant transition from 
the ground state. 
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Figure 1. Nonequilibrium Criterion Curves. The curve for each 
specie-transition represents the "limiting electron 
concentration" for which population fractions of the upper 
electronic state of the transition is 1/2 the equilibrium 
value. Conditions above a curve approach equilibrium, 
while conditions below a curve indicate nonequilibrium. 
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In interpreting Fig. I electron concentrations well above the (N g ) curve 

assure equilibrium population fractions while electron concentrations 

well below the (N g ) indicate collision limited nonequilibrium 

15 -3 

populations. Thus, for = 10 cm the (2+) and the (SR) are in 
nonequilibrium while the (1+) and the (LBH) are in equilibrium. 

In Table 1, a comparison of the population fractions predicted by Eq . 
(5) is given for ^(1+) over a range of temperatures and at two electron 
concentrations . 


TABLE 1 


Population Fractio Prediction of Equation (5) 
for the Upper STa?e of the ^(1+) Transition 


Tg (ev) (N 2 /N x ) 


eq 


°12 N e /A 21 


n 2 /n 


1 


15 13 15 13 

N = 10 N = 10 N = 10 N = 10 

e e e e 

) 


.4 

2.1 

X 

H 4 

o 

1 

OO 

2.8 

X 

l 

o 

2.8 

X 

o 

1 

VJ 3 

1.9 

X 

10 8 

2.5 

X 

io ' 9 

.5 

8.3 

X 

10- 7 

1.4 

X 

10“ 5 

1.4 

X 

10 7 

7.8 

X 

io “ 7 

1.2 

X 

io “ 7 

1.0 

1.3 

X 

10- 3 

3.5 

X 

10 “ 2 

3.5 

X 

10“ 4 

1.2 

X 

io “ 3 

2.8 

X 

io “ 4 

2.0 

5.1 

X 

10- 2 

1.9 



1.9 

X 

10“ 2 

4.9 

X 

10“ 2 

1.4 

X 

io -2 

3.0 

1.7 

X 

io _1 

6.8 



6.8 

X 

io “ 2 

1.7 

X 

10- 1 

4.9 

X 

io “ 2 


In the second column above the equilibrium predictions of Eq. (7) are 
listed. Recall these are also the ratios of the collisional excitation 
and quenching. The third and fourth column list the collision limited 
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solution, Eq. (8), for the two electron concentrations considered. Note 

the ratio of the third or fourth columns to column two is the ratio of 

collision to radiative quenching which appears in the denominator of Eq. 

15 -3 

(5) . An electron concentration of 10 cm is well above the (N g )* limit 

for N 2 (l+) and we find the population fraction from Eq. (5) in column 5 

very close to the equilibrium value in column 2. By contrast, an 

13 -3 

electron concentration of 10 cm is below (N )* for the N.(l+) and the 

e Z 

predictions of Eq. (5) in column 6 are close to the collision limited 
solutions of Eq. (8) in column 4. 

Notice that the collision limited population fraction "prediction" is 
not always below the equilibrium fraction; however, only those values 
below the equilibrium value are valid. Collision limited solutions are 
valid only for N g < (N g )*. Thus, if we make a computation of equilibrium 
and collision limited fractions, we may be assumed that the smaller 
fraction is the valid fraction and Eq. (5) would predict a fraction which 
is below the valid limit. Note that stimulated absorption and cascading 
from higher excited states could raise the fraction above the prediction 
of Eq. (5). 

Now, using the optically thin form of the radiative transfer 
equation, consider the implications of Eq . (5). Combining the two 
equations yields : 

Q/v = hi/A 21 N 2 

■ = h«/(0 12 N e N 1 )/[l+(n 21 N e /A 21 )] (10) 

For the equilibrium limit achieved with N g » (N g )* 
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(Q/ V ) e q = h^A 21 N 2 


(ID 


and for the collision limited with N « (N )* 

e e 

(Q/V) cl = hl/fl 12 N l N e (12) 

The first expression applicable to equilibrium conditions indicates 

the emission from a unit volume is proportional to the number of emitters 

and is limited by the radiative decay rate. The collision limit 

indicates the emission from a unit volume is proportional to the 

excitation rate. Note also that the equilibrium emission scales linearly 

with density of the specie while the collision limit scales as the square 

of the density. When predicting Q/v, it is possible for the equilibrium 

prediction to exceed the collision limiting prediction as it is also 

possible for the reverse to occur. The valid Q/v is the lower of the two 

predictions. The cross over point for Q/v predictions is (N g )*. 

As an illustration of the variation of Q/v with N g , consider Table 2. 

The first column indicates the band system. The electron temperature 

15 -3 

was fixed at lev and the specie ground state concentration was 10 cm 

Column 2 indicates the equilibrium solution. Columns 3, 4, and 5 

indicate the solution of Eq. (10). For the lowest N g only the (LBH) band 

of nitrogen is close to equilibrium with the other transitions well below 

15 -3 

equilibrium. The = 10 cm indicates the (1+) joins the (LBH) in 

equilibrium. At the highest electron concentration considered N g = 10^ 
-3 

cm all states are approaching equilibrium emission. Note the relation 

of these three N 's to the (N )* for the bands of Figure (1). 
e e 

In the above, we have used the optically thin assumption to 
effectively de-couple the population prediction from the radiative 


20 



transfer process while retaining the nonequi librium character of the 
problem. In doing this, we may be under predicting state populations 
slightly. 


Table 2 

Prediction of Volumetric Emission at T = lev and N, = 10^ 3 cm 

e 1 

by Equation (10) for Different Electron Concentrations 

_3 

Transition Q/v watts/cm 


Equilibrium 


N = 10 13 N = 10 15 N = 10 17 

e e e 


N 2 (l+) 

0.39 x 

CM 

1 

O 
» — 1 

0.08 

X 

N 2 (2+) 

0.37 

o 

0.11 

X 

n 2 (lbh) 

0.31 x 

10 2 

0.21 

X 

0 2 (SR) 

260 


0.02 



10 2 

0.37 x 10~ 2 

0.39 

x 10 2 

10“ 2 

0.08 

0.37 


IQ' 2 

0.3 x 10~ 2 

0.31 

x 10~ 2 


2.8 

130 



Next, we shall consider how this simple model can be transferred into 
a simple expression for electronic state temperature while considering 
the vibrational and rotational states of a molecular sources. 
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POPULATION FRACTIONS FOR ELECTRONIC STATES 


Applying the optically thin assumption and considering electron 
impact kinetics to dominate heavy particle kinetics, the expression 
determining the population, N , of the upper state, u, from the rate 
processes between states u and lower state 2 , of a molecule is: 


DN 
u 

Dt 


(Cl„ N.N - A .N N - A .N ) 
0u 2 e u.2 u e u2 u 


(13) 


where Cl. and Cl . are the kinetic rate coefficients for excitation and 
2u u 2 

de-excitation of state u by electron impact and A ^ is the spontaneous 
emission coefficient between u and 2. 

The total population N^ of the upper electronic state which is 
designated by "2" is obtained by summing over all the upper states 
associated with the upper electronic level. So 


DN„ DN 
— ^ = S — H = N 
Dt u Dt 


ZEfl. N 
. 2u e 
u 2 


-N SSO.N-SSA.N 
e . ul u . ut u 
u 2 u 2 


(14) 


Considering that the upper state u is designated by the electronic state 
number 2 and the vibrational state number v' and that the lower state 2 
is designated by the electronic state number 1 and the vibrational state 
number v' ' , it is clear that the summations in Eq. (14) are over all the 
v' ' and the v' associate with the (2,1) transition. Although the 
rotational structure of the states has not been considered, to be 
consistent with the bandless model used with the radiative transfer 
equation, the summation over rotational states would yield expressions 
identical to those below. 
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Now using the Boltzmann distribution for the fraction f of the 
population of an electronic state which is in the specific vibrational 
state v' 


V N 2 = f V = »2V /N 2 = ex P < - i v' /kT v2 )/Q 2 <T v2> <l5a) 

and f , , for the lower electronic state fraction 
v' 

= f v , , = = ex P (-% (> AT v i)/Qi(T v i) (15b) 

Where the Q's are the vibrational partition functions for the respective 
electronic states which are evaluated at the Boltzmann vibrational 
temperatures: for the upper electronic state and for the lower 

electronic state. 

Recall that the spontaneous emission coefficient can be factored into 
an electronic transition probability A ^ and a Franck-Condon factor 

Vv' ' so that 


A ul = A 21 V 


(16) 


Thus, using Eqs . (15) and (16), terms such as 


2 2 A „N = A 01 N_ 22 q , , , f , 
. u! u 212 , , , v’v' ' v' 

u 2 v' v' 


(17) 


The sum of the Franck-Condon factors over any column or row is unity so 
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The sum of the vibrational state population fractions over all upper 
electronic state vibrational levels is also unity so the double sum on 
the right side of Eq. (17) is unity and the equation reduces to: 


2 2 A u* N u = A 21 N 2 
u e 


(18) 


Similar arguments can be used for the collisional rate coefficients so 
Eq. (14) can be re-written as: 


DN 

— - = o 
Dt 


12 N l N e 


n 21 N 2 N e A 21 N 2 


(19) 


Where and are the populations of the lower and upper electronic 
states and and A 21 are t '* ie c °ll^ s ^ ona l excitation rate 

coefficient, collisional quenching rate coefficient, and the radiative 
(spontaneous) quenching rate coefficient. 

The steady state population ratio is achieved by setting Eq. (19) 

/ 

equal to zero so: 


V N 1 = n !2 N e /(A 21 + n 21 N e ) (20) 

which is identical in form to the illustration in the previous section 
The equation will be written as: 


n 2 /n x = (R 1 R 2 )/(R 2 + r l ) 


( 21 ) 


where r = A 01 ^ is the radiative lifetime and R, = rfl, „N and 
2 1 I 12 e 
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R 2 = <VVeq “ <*2V T v2>'*lV T vl» "P <“ 12 '“. > = n i2 /n 2l 

The quantity R is ratio of collisional excitation to radiative quenching 
rates, the quantity R£ is the ratio of upper to lower electronic state 
populations in equilibrium*, and the ratio *- s t ' ie rat i° of 

collisional to radiative quenching rates. 

Rather than an expression for populations such as Eq. (21), it is 
more useful in radiative transfer applications to express a population as 
an electronic state temperature, 

T E2 = <e 12 /k) ln t (N l/ N 2 ) <W (6 2 /g l )] (22) 

Although developed for a resonant transition, i.e. 1 is the ground 
state and 2 is an excited state, radiative quenching to other excited 
states may be considered simply by substituting an appropriate radiative 
lifetime for r in Eq. (21). 

In closing, note the relationships between individual state rate 
constants and electronic state rate constants are: 


^u 2 ^21^v' v' ' ^v' 


A ufi q v'v' ' f v ,//r 


(23) 


* Note equilibrium is in terms of the electron temperature T and the 


upper and lower electronic states vibrational temperatures and 


T vr 
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RADIATIVE TRANSFER EQUATION 


This section is devoted to an examination of the radiative transfer 
equation and shall show how the standard or classic form of this equation 
in terms of source function and absorption coefficient can be constructed 
for nonequilibrium conditions from a modification of the conventional 
absorption coefficients and source functions computed by equilibrium 
computer codes. Solution procedures for the classic radiative transfer 
equation will not be considered as adequate treatment is available in the 
form of computer codes and their documentation. Similarly, we will not. 
review the computation of absorption coefficients as the procedure for 
computing the bandless model is incorporated into the same codes. We 
shall concentrate on the necessary modifications and their logical 
development. The development is easily extended to more complex models 
such as the JORL. 

Assuming that either excited state population fractions or 
temperatures are available from solutions of the state population rate 
equations, the starting point is the "microscopic" radiative transfer 
equation expressed in terms of the Einstein relations for bound-bound 
transitions. The equation, Ref. [18], for the rate of change of 
radiative intensity I in the frequency range Ai/ at frequency u, in the 
direction z per unit solid angle, considering the processes of 
spontaneous emission, stimulated absorption and stimulated emission 
between an upper electronic state 2 with population ^ and a lower state 
1 with a population is: 


M — = <hi//47r) S (A^i ) N^ l) ^ l) - I J ) N 1 (l) * 1 (l) )] (24) 

dz i 


>12 "1 


21 2 


26 




where the summation over i includes all species and bands which 
contribute to the Au interval and where and w bich may be though of 
as "line" shape functions for the transition, are: 

(25a) 
(25b) 

where the sums are over all vibrational-electronic transitions with 
centers in the interval Au . Recall the quantity fi is the direction 
cosine between the intensity ray and the direction z. The above 
expression is consistent with the bandless model which assumes the 
interval Au to be small enough so that properties are constant in the 
interval but broad enough to encompass all of the rotational lines 
associated with the transition. 

Recognizing that the terms multiplied by 1^ comprise the absorption 
coefficient including stimulated emission: 

// = S jP (1 - «£ x) ) (26a) 


*2 = l Vv'^v'/^ 

Au 


= 2 Vv" f v" /Al/ 
Au 


with 


5 

u 


<i)_ 


< B 21 N 2V B 12 N lV 


(i) 


(26b) 


and the absorption coefficient 


/*„ (l> = (W4 jt) (B 12 N iV (1) 


(26c) 
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Factoring this term in Eq. (24) yields the standard form: 


dl 


dz 


= n (S - I ) 

V V V 


(27a) 


where the source term takes the form: 


S = S (A. 1 N 1 $ 1 ) (l) /^' 

i / . 21 1 1 v 

l 


(27b) 


With the use of the relationship between the Einstein coefficients: 


3 , 2 


A 21 = (2hi/ /c ) B 2] _ ; B 12 = (g ^g p 21 


(28) 


the terms 


3 , 2, 


A 21 N 2*2 = < B i 2 N 2V < 2ht/ > (S 1 /S 2 ) 


(29a) 


which comprise the sum and the numerator of the source function. Now 
multiply and divide by which are evaluated at the lower electronic 

state vibration temperature yielding: 


3 , 2 


(B 12 N lV = < 2lu/ / c > ( Sl N 2V g 2 N lV 


(29b) 


with the introduction of a terra, Eq. (26b), 6^ = 2 ® 2^2^ 1^1’ t * ie a ^ ove 

can be written as : 


[(B 12 N iV (1 “ K )] [< 2 ^V> S u /a - K )] (29c) 
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which one recognizes as the product of the absorption coefficient and 
the Planck or Blackbody source function provided 6 reduces to the simple 
exponential of frequency and temperature . The source function becomes : 


s <K )W ( V (i>/ 'C (30) 

i 

where the absorption coefficient for the individual band interval for the 
i'th transition: 


(p') (i) = (hu/Un) (B 12 N 1 $ 1 ) (l) (l - 5^ (1) ) 


and the source function for the individual band interval from the i'th 
transition is: 



(2hi/ 3 /c 2 ) « (l) /( 1 - 


5 (i) 


) 


with 


= (S 1 N 2 ^ 2 /S 2 N i^i )(i) (3D 

First note that the ratio of upper electronic state population to 
lower electronic state population: 

N 2 /N 1 = (g 2 / g]L ) (Q 2 ( T v 2 ) /V T vl ) exp (- £ 2 / kT E2 )/ ex P (- £ l/ kT Ei> ( 32 > 

which is not a simple function of temperature and frequency. 


29 


Furthermore, re-examining the sums of and $> 2 > Eqs. (25 a & b) , it is 
not obvious that these sums over the Ai/ band interval are equal other 
than under the conditions that Ay encompasses the entire band so that the 
sum encompasses all upper and lower states. To explore the relationship 
between and $ 2 > it i s convenient to recall that for frequencies u in 
the Ay interval are represented by the upper and lower state vibrational 
energy terms satisfying: 


ht/ 


: 12 


+ € 


(33) 


which can be used to solve for e , and used with Eqs . (25) to yield a 
ratio 


*3/*! = 'W/Wl [ V T v2 )/Q 2 <T v2 >)e ' h< ‘''''°o )/kT ''2 ° 4) 


where y = 6-,,,/h. 
oo 12 

Combining Eqs. (34 and 32), we find that 

6 V = K*^) '/<*]_<*!>"] R exp (-hi//kT v2 ) (35a) 

where the double prime indicates the functions are evaluated at and 

the single prime indicates the functions are evaluated at and 


R = exp (~£ 12 /kT t )/exp ( £ 12 /kT v2 ) 


where T , 


the transition temperature is 


(35b) 
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(35c) 


T t T E2' / ^ 1 + (€ l //e 12 ) (T E1 T E2 ) / T E1^ 

Note first that for equilibrium R = 1 and the partition function and 
line shape function ratio is unity so we find S = exp (-e^/kT), the 
familiar Planck exponential. Next, note that in general, for 
nonequilibrium S is a function of 4 temperatures, i.e. a vibrational and 
electronic state temperature for both the upper and lower states of the 
radiative transition. If the transition is to a ground state, = 0 and 
the transition temperature reduces to the upper state electronic 
temperature; however, both upper and lower vibrational state temperatures 
are required. Next, note if all upper and lower states share a common 
value of vibrational state temperature then the line shape and partition 
function ratio terms reduce to unity. Clearly, under most nonequilibrium 
conditions , the S ^ term does not reduce to the Planck exponential ; 
however, it can be re-written as a frequency shifted Planck term 
evaluated at the upper state vibrational temperature. 

S^ (l) = exp (~h(i/ + A (L) )/kT v2 (l) ) (36a) 

where the extent of the frequency shift is : 

A (l) = *' 00 [(T v2 -T t )/T (; ] + kT v2 In [ (Q^) " /(Q^) ' ] (36b) 

Thus, the algorithm for the standard form for the radiative transfer 
equation, Eq. (27a) consisting of the absorption coefficient Eq. (26 a 
and c) and the source function by Eq. (30) with & given by Eq . (36). 


31 



In adapting this algorithm to existing programs, the absorption 

coefficient could be taken care of by an equilibrium code provided the 

T , is used in its evaluation. Note T_, is also necessary if the lower 
vl El 

level is not a ground electronic state. The 6 would have to be 
evaluated by using the upper vibrational state temperature and computing 
the frequency shift which is a function of all four temperatures 
associated with the transition. 

Note that the shift in the Planck exponential , terms are positive for 

T „ > T - and T „ > T . This shift would attenuate UV and increase the 
v2 vl v2 t 

potential contributions in the IR. High electronic state temperatures 
would reverse this pattern. 

The reduction of the above to the optically thin approximate, such as 
discussed in earlier sections, would indicate the dependence upon T^ and 
T is not as strong as the above expressions seem to indicate. 
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CONCLUSION 


The objective of this study has been to develop the algorithm which 
will allow nonequilibrium flowfield codes to be coupled to radiative 
transfer codes which are currently limited to equilibrium applications. 
The algorithm consists of Eq. (22) for the prediction of electronic state 
population fractions (electronic state temperature) and Eqs . (35c), (36), 
(26c) , (26a) and (30) for the compution of spectral absorption 
coefficients and source functions for use in Eq. (27a). These 
expressions represent a thorough treatment in the context of the bandless 
model which in the past has proven adequate for equilibrium radiative 
transfer predictions. Perhaps the approach is more than adequate as the 
current algorithm requires four temperatures for each transition. 

Clearly, a more sophisticated model can be developed; however, it would 
require more information ie. rotational temperatures of the electronic 
levels . 
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